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Comprehensive Thermodynamic Theory of Stability 
of Irreversible Processes (CTTSIP): The Setup for 
Autonomous Systems and Application 

Vijay M. Tangde*, Seema G. Rawat and Anil A. Bhalekar f 


Abstract 5 In the present communication, for the first time 
a version of Bhalekar’s Comprehensive Tl dy z Theory 
of Stability of Irreversible Processes (CTTSIP) namely CTTSIP 
for Autonomous Systems has been described. The Autonomous 
Systems are defined in terms of the differential equations of 
motion in which the time variable t does not appear explicitly 
and hence all functions remain only implicitly time dependent. 
As in CTTSIP one uses the fabric of Lyapunov theory of 
stability of motion, there evolve the conditions of 
thermodynamic stability, thermodynamic asymptotic stability, 
thermodynamic stability under constantly acting small 
disturbances and thermodynamic instability. As an example 
this new version of CTTSIP has been applied to chemically 
oscillating system of Brusselator Model. The regions of 
thermodynamic asymptotic stability and stability under 
constantly acting perturbation emerge. 

Index Terms —Thermodynamic Stability, Autonomous 
Systems, Irreversible Thermodynamics, Chemical Oscillations, 
Brusselator Model. 

L Introduction 

One of the present authors (AAB) has developed a 
Comprehensive Thermodynamic Theory of Stability of 
Irreversible Processes (CTTSIP) [1-4], The basic fabric of 
CTTSIP is that, it is woven using Lyapunova second (direct) 
method of stability of motion [5-8] and the second law of 
thermodynamics [9-12]. Indeed, CTTSIP is compact and 
similar to the Gibbs-Duhem stability theory [9-11, 13] of 
states of equilibrium thermodynamics. Unlike the Glansdorff 
and Prigogine theory [10] for so-called local equilibrium 
states, in CTTSIP, 

1. one does not require the assumption of local 
equilibrium, (in fact, it has been demonstrated 
elsewhere [1-4] that its is a physically non-existent 
assumption), 

2. importantly it has a sound thermodynamics basis [1, 
2, 4], 

During the last decade, this framework has been well tested 
and applied to a variety of real irreversible processes from 
laboratory [1, 4, 14-18] to industry [19, 20]. Now on taking a 
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stock of the status of CTTSIP we realized that the majority of 
thermodynamic systems fall under the category of 
Autonomous Systems as described by the original stalwarts 
who propagated Lyapunova second method of stability of 
motion (see for example: [5, 21]). Hence, we take an 
opportunity to describe, in this paper, a version of CTTSIP 
applicable to Autonomous Systems and illustrate its utility by 
applying it to a representative autonomous system, namely, 
oscillatory chemical reaction of Brusselator model. In 
Appendix A we have briefly described Lyapunova direct 
method of stability of motion for ready reference and the 
original setup of generalized CTTSIP has been presented in 
Appendix B. 

IL CTTSIP Setup for Autonomous Systems 
There are four building blocks of CTTSIP setup, 

namely 

1. the identification of thermodynamic coordinate and 
corresponding space, 

2. defining of thermodynamic perturbation coordinates, 

3. construction of thermodynamic constitutive equations 
in terms of perturbation coordinates and 

4. defining of thermodynamic Lyapunov function that we 
describe in the following subsections. 

A. Thermodynamic Perturbation Coordinates, Space and 
Constitutive Equations 

Let X? (t) a be the thermodynamic coordinates of 
the given process on the unperturbed trajectory and X t (t) a 
be that on the perturbed one. The thermodynamic 
perturbation coordinates , OC i (t) a then get defined as, 
a i (t) = (x i (t)-x?(t))^Q, (i = 1 , 2 , •••,«) (2.1) 

and hence the equation of unperturbed trajectory is obtained 

a? (f) = 0 (i = 1,2,(2.2) 

where in the above equations the superscript 6° Odenotes the 
quantity on the unperturbed trajectory. Further, as per 
requirement of Lyapunov theory the perturbation coordinates 
should satisfy the following condition namely, 

I a t (0)1=1 X,. (0) - x° (0) \=\ (i = 1,2, •••,«), (2.3) 

Where, is a sufficiently small number 
where for the sake of brevity we have shown by O 0 6 the 
values of the parameters at t = t 0 , the time at which the 
perturbation has been effected. 
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In this way Gift) a are the sufficiently small 
perturbation coordinates in the domain 
t>t 0 , t> 0, \a i (f)\<£, tf>0, (2.4) 

where (f is a sufficiently small positive number. 

These perturbation coordinates in turn determine the 
thermodynamic perturbation space. 

The thermodynamic perturbation constitutive 
equations of motion within the perturbation space in general, 
are described as, 

< ^ L = f(a 1 (0, a 2 (0, « 3 (0, • • •, a n (t)). (2.5) 

at 

Thus, we see that the trivial solution of equation (2.5) is 
equation (2.2). 

B. The Choice of Thermodynamic Lyapunov Function 

For defining thermodynamic Lyapunov function we 
take the help of Clausius-Duhem inequality [12, 22-24], 
namely 
ds 

p -1- V • J s - a s > 0, (2.6) 

dt 

where p is the mass density, S is the per unit mass 
entropy, J s is the non-convective entropy flux density and 
<7 S is the entropy source strength. 

In CTTSIP[l-4] we use the above described sign 
definite entropy source strength function and define 
thermodynamic Lyapunov function, 4 , as, 
4=(<7 s (0-<X°(0)^0, (2.7) 

where, cr (t) and <7°(f) are the respective entropy source 
strengths and are positive definite quantities[12, 25, 26] as 
guaranteed by the second law of thermodynamics appears in 
the Clausius-Duhem inequality[12, 22-24], Since we are 
using the thermodynamic function, namely the entropy 
source strength in defining the Lyapunov function it 
automatically gets thermodynamic sanction. Hence, C s is 
the excess rate of entropy production per unit volume. 

Thus the self consistency of this choice of 
thermodynamic Lyapunov function gets demonstrated by the 
fact that when it is the case of thermodynamic stability of an 
equilibrium state the second law of thermodynamic 
guarantees <7° = 0 and hence the thermodynamic 
Lyapunov function reads as 4 = (T s (t) , that is the entropy 
source strength itself takes the role of Lyapunov function in 
this case. Which is nothing else but the Gibbs-Duhem 
stability theory of equilibrium states [10, 27]. 

Recall that, the thermodynamics of irreversible 
processes provide the following dependence of entropy 


source strength [12, 25, 26], <7,(0 , namely, 

cr (t) = <7 s (Xi,X 2 ,X 3 ,. ,X n ) (2.8) 

or more elaborately as, 

cr (0 = cr.(x, (t),x 2 (t),x 3 (t), . ,x n (t)) (2.9) 

for autonomous systems. 


Therefore, the fundamental dependence of C s is 
obtained as, 


4 (0 = C s (a, {t), a 2 (0, <Z 3 (0, -,a n (/)), 

or more elaborately as, 

4(0 = C,faft),a 2 (t),a 3 (t),-,a n (t))> 
e(a v a 2 ,a v --) > 0 for t>t 0 


4 (0 = 4 («! (0. a 2 (0, «3 (Or • •, a n it)) 

( 2 . 12 ) 

<-Tj{a v a 2 ,a 2 ,---) < 0 for t>t 0 , 
corresponding to the two options of equation (2.7), where 
4 is differentiable function [28] and £ and Tj are 
continuous positive numbers. Thus, from equation (2.10) we 
obtain the following equation of unperturbed trajectory, 
namely: 

4(0 = 4 (0,0,0,. ,0)e0 for t>t 0 , (2.13) 


in both the cases of equations (2.11) and (2.12) with 

<£•(0,0,0,.) = 0 and also 

(2.14) 

77(0,0,0,.) = 0. 

From the above description it gets established that the 
majority of thermodynamic systems constitute the 
Autonomous Systems as has been defined in Lyapunov theory 
of stability of motion [6, 8], 

The total time derivative of thermodynamic 
Lyapunov function, 4(0 > i n the case of Autonomous 


systems read as, 

dC,_y( 541 da, 
dt ~y\da i ) dt 


(2.15) 


C. Thermodynamic Stability of Autonomous Systems 


1. For autonomous systems along with the condition 
of equations (2.11) and (2.13) if, the following condition is 
followed, namely: 



then unperturbed (real) trajectory under investigation is 
guaranteed as thermodynamically stable. 

2. Whereas along with the conditions of equations 
(2.12) and (2.13) if, we have. 


4(0= 


dC s 

dt 



(2.17) 


then also the unperturbed(real) trajectory under 
investigation is guaranteed as thermodynamically stable. 

Recall that the fundamental dependencies of entropy 


source strength on OC^S read 


cr = cr (a x ( t ), a 2 (t), a 3 (t), .) > 0, (2.18) 

cr° = cr,(0,0,0,0,.) > 0 (2.19) 

and the positive definiteness of (2.18) and (2.19) is 
guaranteed by the second law of thermodynamics [10-12], 


D. Thermodynamic Asymptotic Stability 

On following the Lyapunov theory we arrive at the 
following description of thermodynamic asymptotic stability 
of a real trajectory in the said two cases namely; 

1. The validity either of equations (2.11) and (2.13) 
and in addition [21] 
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Jjxf>p 2 > 0, t>T 0 >t 0 

and 


4 ( 0 = 


^4 

dt 



< -p < 0. 


( 2 . 20 ) 


( 2 . 21 ) 


2. Equations (2.12), (2.13), (2.20) and 


4 ( 0 = 


dC s 

dt 



>P> 0, 


( 2 . 22 ) 


with d(X i / dt given by equation (2.5) in the domain 
prescribed in equation (2.4) where 

P = P(a x ,a 2 ,a 3 , . ,a n )> 0, (2.23) 

4° =4(0,0,0,.,0) = 0, (2.24) 

that is, P is a strictly positive number and vanishes only at 
the origin or equivalently on the real trajectory. This in terms 
of 4 means that d£, s / dt vanishes only at the origin or 
equivalently on the real trajectory and outside it is strictly 
negative (positive) quantity. 


E. Thermodynamic Stability Under the Constantly Acting 
Small Disturbances 

In addition to the conditions prescribed in either 
equations (2.11), (2.13), (2.20), (2.21), (2.23) and (2.24) or 
equations (2.12), (2.13), (2.20), and (2.22)-(2.24) if the 
partial derivatives djO, / 5(2; remain finite, then 
thermodynamic stability guaranteed, is also against the 
constantly acting small disturbances. This assertion is based 
on MalkinQ theorem[6, 8, 21], 


III. Theorems of Thermodynamic Stability for Real 
Trajectories of Autonomous Systems 
The above presented thermodynamic stability 
description for Autonomous Systems culminates into the 
following theorems, namely: 

Theorem I. For a system of equation (2.4) of the 
perturbed motion in the thermodynamic perturbation space 
determined by equations (2.1) - (2.3) and restricted by 
equation (2.5), if there exists a differentiable thermodynamic 

Lyapunov function 4(0 = 4 {^a x (t),a 2 (t),a 3 (0>.) 

defined by equation (2.7), which satisfies the following 
conditions in the neighborhood of the coordinate origin for 
t > t 0 ; namely: 

4(0 = £, (Oj(0,« 2 (4« 3 (0,. ,a n (t))> 

/ 3 C«) 

s[a x ,a 2 ,a 3 , . ,cc n )> 0, 

with 4 (0 = 4 (0,0,0,.,0) = 0 , where E is a 

strictly positive continuous function that vanishes only at the 

origin, E° =<£■(0,0,0, . ,0) = 0, that is, 4(0 has a 

strict minimum at the origin, and the derivative 



then the trivial solution of the system of equation (2.4), 
namely (2 ; ° = 0, that is equation (2.2), constitutes a stable 
motion in thermodynamic space. Similarly, the pair of 
equations, 

4(0 = 4 (aft),a 2 (t),a 3 (t\ . ,a n (t))< 

-rj{a x ,a 2 ,a 3 , .,«„)< 0 





for t>t 0 establish that the trivial solution of the system of 
equation (2.4), namely OC i = 0 , that is equation(2.2), 
constitutes a stable motion in thermodynamic space. 

Theorem II. For the system of equations (2.4) of the 
perturbed motion in the thermodynamic space defined by 
equations (2. l)-(2.3) in the domain determined by equation 
(2.5), if there exists a differentiable thermodynamic 
Lyapunov function, 

4 (0 = 4 { a i ( 0 » a 2 ( 0 » a 3 ( 0 > . > a „ (0), defined by 

equation (2.7), which satisfies the following conditions in the 
neighborhood of the coordinate origin for t>t 0 ; namely: 


4(0 = 4 ( a n oc 2 ,a 3 , ., a n ) > 

/ 3 tfW) 

s[a x , a 2 ,a 3 , . ,a n )>0, 

with 4 (0 = 4 (0,0,0,., 0) = 0 , where £ 

is a strictly positive continuous function that vanishes only at 
the origin, £° =<£'(0,0,0,. ,0^ = 0, that is, 4(0 has 


a strict minimum at the origin, and the derivative 


4 ( 0 = 


44 

dt 



(H-2) 


where £ 2 is strictly positive continuous function that 
vanishes only at the origin, E 2 — E 2 (0, 0 , 0 , • • ••, 0) =0 , 
then the trivial solution of the system of equations (2.4), 
namely (X- = 0 , that is equation (2.2), constitutes an 
asymptotically stable motion in thermodynamic space. 

Similarly the following pair of equations, namely: 

4(0 = 4 ( a va 2 ,a 2 , ., a n ) < 

-r] (cij, a 2 , a 3 , ., a n ) < 0, 


4(0 = 


44 

dt 



da j 

'~di 


> 7 2 > 0 , 


for t>t Q establish that the trivial solution of the system of 
equation (2.7), namely Ot l = 0 , that is equation (2.2), 
constitutes an asymptotically stable motion in 
thermodynamic space. (II-3) in addition if the derivatives 

54 /5(2; (z = 1,2,3, . ,n) are finite (bounded in 

absolute value), then the said trajectory (real) in 
thermodynamic space is obtained as stable under the 
constantly acting small disturbances. 
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IV. Application to chemical oscillations in 
Brusselator model 


We exemplify the use of CTTSIP for autonomous 
systems described in sections 2 and 3 by applying it to the 
chemical oscillations in Brusselator model. In the 
mathematical model of Brusselator two intermediate species 
show oscillating behavior. The schematic representation of it 
is [29, 30], 
h 

A^>X (4.1) 

2X + 74-3X (4.2) 


B + X 


■Y + D 


(4.3) 


xXe 

where A and B 


(4.4) 

the reactants, X and Y are the 


intermediates, D and E are the products. The rate equation 
for such autocatalytic reactions are fundamentally non-linear. 



Figure 1: Schematic representation of Brusselator model 

In the Brusselator model, it is assumed that the 
concentration of reactants cAOand d30are kept constant at a 
desired feeding rate by means of continuous-flow stirred tank 
reactor (CSTR) as illustrated in the figure 1 and the 
concentration of CDOand CEOare removed continuously from 

k k k k 

the system as they are being produced. 1 , 2 , 3 and 4 

are rate constants of the respective chemical reactions. 

x/mol dm -3 & y/mol dm -3 



Figure 2: Oscillating behaviour of d<0 and d/0 with time, 
where the bold line indicate oscillations in the concentration 
of intermediated-^ Oand other line represent the oscillations 
in concentration of CAOintcrmcdiatc. 

From chemical kinetics [31], the rates of change of 
intermediate species on unperturbed trajectory are given by. 


— = k.A 0 + k 2 {X°f Y° - k 3 X°B° - k A X°, 

dt (4.5) 

— = - k 2 (X° f Y° + k s X°B°, 

dt (4.6) 

and that on the perturbed trajectory by, 

—=k.A° + UX) 2 Y-LXB° -LX, 

dt (4.7) 

-=-k 2 (X) 2 Y + k 2 XB°. 

dt (4.8) 

Notice that as the concentrations of and are 
maintained constant in the reactor by the constant rate of feed 
in, the concentrations of A and have been used in 
equations 4.5 to 4.8. Also, notice that for the sake of bravity 
we have shown the concentrations of chemical species by 
o 

their chemical symbols itself. The superscript on the 
symbols describe the concentrations on the unperturbed 
trajectory and without superscript they represent 
concentrations on the perturbed trajectory. In this case, 
oscillations [13, 29, 30] sustain as long as reactants are 
continuously fed into and products are removed from the 
reactor. The oscillating behavior of intermediates 6 X Oand 
&iO are generated graphically by subjecting their rate 
expressions for numerical computations using Mathematica 
Software and is represented in figure 2. 

Let us now consider a case when concentrations of both 
autocatalytic intermediates that oscillate, namely, X and Y 
are perturbed simultaneously by sufficiently small amount 
say afiX oand , that is, 

SX=X-X° and JY = Y-Y°. (4.9) 

The equations of motion within the thermodynamic 
perturbation space are then obtained as, 
d(X-X°) _ d(SX) 

dt ~ dt (4.10) 


= ( k 2 Y° (2X - SX) - Bk 2 -kf)SX 

and 

dt dt (4.11) 

In the present case, the operative expression of 

£ 

thermodynamic Lyapunov function, 5 , would be of the 






15(^7 )) 


(SY)^O 


and the operative expression of total time derivative of 
C , would be 


(4.13) 

s of entropy 



d(SX) 

f JO 

d(SY) 

[d(SX)) 

dt + 

[d(SY)) 

dt 


Let us recall the definition of global level r; 
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production [9, 11, 32-33] on the unperturbed trajectory as, 

z o = A° djfl { A? d(j)l , 

T dt T dt (4.14) 

d({t , X ^ Q 

T dt T dt 

and accordingly on perturbed trajectory it would reads as, 

s = A <*(£),- , A <*(£)„■ | 

T dt T dt (4.15) 

A,. ( A, <*(#)*, Q 

T dt T dt 

On substituting the rate expressions and further 
solving, we get, 

_A 

T 

A A 

(LXB°) + YXl (LX) > 0. 

T 3 r 4 

Notice that the positive signs of equations (4.14), (4.15) and 
(4.16) are dictated by the second law of thermodynamics. 
Using CTTSIP [4] for Autonomous Systems, the 

thermodynamic Lyapunov function, C s , is constructed by 
differentiating equation (4.16) partially with respect to ^X 
first then with respect to ^Y , we get, 

T 


= A A 0 )+A (k 2 X 2 Y) + 


0(JX) 




(4.17) 


7? \L 


LXB°. 


d(SY) {T 

On imposing the condition of real trajectory, it gives 


f-^-1 =-^\A° +[2^-R\k 2 X°Y 0 + 
J X° 1 { T J 2 

= ^+^ 2 (. 




f —1 

l d(SY)J 


f 2 (^°) 2 - “ \k 3 X°B°. 


(4.20) 


Accordingly, the operative expression of thermodynamic 

4 


Lyapunov function, 

4=1 


/ as, \ 

V )+ f.®4 

l d(SX)J 



would be of the form, 

(^)^O 


on substituting equation (4.17) and (4.18) in the equation 
(4.21), expression of thermodynamic Lyapunov function, 
C s , gets simplified as. 


[)Y +R ) k A+[)y+ R Y] (SX) 


Now, the operative expression of total time derivative of 
C s , would be 

dC s J 8C S \ d(SX)J a4 \ d(SY) 

M v d(SY) J' 


dt 


d(SX) 


dt 


dt 


(4.23) 

On partial differentiation of (4.22) with respect to the 
perturbation coordinates SX and SY , we get the 
expressions of gradient of C s as, 

d(SX) X° " 1 " \ T 
^+ R j k 3 B° + f ^ +R j k A is finite 


-^={jr +R ^k 2 ( x y -{ R ) k . x ° B0 ™ finite. 

(4.25) 

On substituting equations (4.10), (4.11), (4.24) and (4.25), 
in equation (4.23) we get, 



(4.26) 


A. Thermodynamic Stability Discussion 
Because For determining the thermodynamic stability of 
Brusselator model the final expressions of C s and df / dt 
given in equations (4.22) and (4.26) are used. Now, in view 
of the complex nature of both the expressions, they are solved 
computationally using software Mathematica 5.0, developed 
by WOLFRAM RESEARCH INC., USA. The schematic 
presentations of this model have been generated using this 
software. 

Following numerical constants have been used for 
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computational inputs, namely, k x = lmol dm' 3 s' 1 ; k 2 = lmol 
dm' 3 s' 1 ; k 3 = I mol dm' 3 s' 1 ; k 4 = lmol dm' 3 s' 1 ; A =1 mol 
dm' 3 ; B =3 mol dm' 3 ; Y =1 mol dm' 3 and X = 1 mol dm' 3 . 
Results of computations are depicted graphically in figure 3 
and figure 4. Figure 3 shows smooth and continuous decrease 
of positive C s while the figure 4 of dC s / dt shows its 
negative values which are not continuous and this may be 
attributed due to feedback mechanism of autocatalytic 
species leading to oscillations. This means that, natures of 
both the graphs favor the thermodynamic stability. Now, as 
per the description in CTTSIP above in section 2.4, the 
thermodynamic asymptotic stability of any process is 
guaranteed if in addition to opposite signs of C, s and 
dJZ s / dt , the following condition remains satisfied, namely, 

Y±. SX fi ^Px>° and ~ Py > °> 

~ (4.27) 

where, p x , p Y are any arbitrary positive definite 
constants. Further, to establish that the process is 
thermodynamic stable under constantly acting small 
disturbances as per Malkin’s theorem, we have verified the 
finiteness of the gradients of >C s with respect to its 
perturbation coordinates computationally. Corresponding 
graphical variation of gradients have been shown in figures 5 
and 6. Figure 5 shows that the gradient of thermodynamic 
Lyapunov function C s , with respect to perturbation of 
intermediate &0 initially, decreases a while and then starts 
increasing with progress of time, reaches certain maxima and 
further decreases back to original level. Moreover, during the 
course of this trend of decrease, increase and again decrease, 
the gradient remains finite and there is no sign of its 
deflection towards infinity. Similar behavior is observed in 
the case of gradient of thermodynamic Lyapunov function 
£, s , with respect to the perturbation coordinate of 
intermediate disestablishing the finiteness of the gradients. 
Thus, for Brusselator model, the thermodynamic stability 
under constantly acting small disturbances is also 
established as per Malkins theorem [6, 8], 



Figure 3: Plot of C s with time when the concentrations of 
intermediates are perturbed. 


cL£,/dt/Jd* 3 k-V ; 



Figure 4: Behavior of d£, s / dt with time when the 
concentrations of intermediates are perturbed. 

3Z x d:SX; 



Figure 5: A plot of dC s / d(SX) with time showing 
ftniteness of C, s against the perturbation in X. Therefore, the 
process is thermodynamic stable under constantly acting 
small disturbances as per Malkin's theorem [6, 8]. 

&£ J ,/a(3¥)/JK‘s-Wl J 



Figure 5: A plot of dC s /d(SY) with time showing 
ftniteness of C, s against the perturbation in Y. Therefore, the 
process is thermodynamic stable under constantly acting 
small disturbances as per Malkin's theorem [6, 8], 


V. Conclusion 

In the present paper, we have reported for the first 
time the version of Bhalekara Comprehensive 
Thermodynamic Theory of Stability of Irreversible Processes 
(CTTSIP) for the Autonomous Systems. Since in autonomous 
systems the explicit time variable does not appear, hence all 
the functions remain time dependent only implicitly. 
Therefore, in the CTTSIP for Autonomous Systems the local 
time derivative of thermodynamic Lyapunov function does 
not appear. Additionally, to guarantee that the process is 
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thermodynamically asymptotically stable one needs to satisfy 
an additional conditions described in equation 2.20, that is a 
crucial ingredient of the theory. In the present 
communication, we have demonstrated the use of the 
thermodynamic Lyapunov function as defined in CTTSIP 
also for Autonomous Systems in a systematic way to predict 
the thermodynamic stability aspects of the processes. As a 
representative example, we have discussed, the 
thermodynamic stability of Brusselator model which consists 
of a set of autocatalytic oscillating reactions. Appropriate 
perturbation coordinates are selected within the 
thermodynamic space namely, SX and SY 
Corresponding thermodynamic Lyapunov function is 
constructed as per the steps of CTTSIP and the expression of 
their total time derivative was obtained. Finally, by 
generating the graphical representations of these expressions 
using the Mathematica Software, the thermodynamic 
asymptotic stability and stability under constantly acting 
small disturbances get demonstrated in the present case. This 
is the reason that the intermediates continue the usual 
oscillations even after perturbation in its observed natural 
trajectory. Evidently, from all our earlier studies using 
generalized CTTSIP and present study using CTTSIP for 
Autonomous systems, it gets lucidly illustrated that there are 
no restrictions at all on the applications of CTTSIP of 
Autonomous Systems. Thus the comprehensive character of 
the framework of CTTSIP of Autonomous Systems get 
adequately illustrated in this paper. 

Appendix 

A. Lyapunov’s Direct (Second) Method of Stability of 
Motion 

For the convenience of readers herein we describe a 
gist of Lyapunova direct (second) method of stability of 
motion. The reader is advised to refer the original references 
on this subject cited in this paper to get its details. 

1. Let the given differential equations of the perturbed 
motion be, 

A-H****..■*■> <»-*•«•. 

(6.1) 

The trivial solution of equation (6.1) is, 

x t = 0 (i = 1,2,3,. ,n ), (6.2) 

where X { a have been defined as, 

*/(0 = lv ; ( 0-^°(01 

x i (0) =1 y t (0) - jv,° (0) I=(6.3) 

where y-(t) are the coordinates of the real motion under 
investigation and correspondingly y i (t) are those for the 
perturbed trajectory. The equation of unperturbed motion 
then reads as, 

x’(t) = 0 and x°(0) = 0. (6.4) 

Thus X t a are the small perturbation coordinates in the 
domain, 

t>t 0 ,t 0 > 0, X t <H, H> 0, (6.5) 

where H is a sufficiently small positive constant and 


X (t,0,0,0, .,0) = 0 (/ = 1,2,3,., n ). (6.6) 


2. Let V(t,X l ,X 2 ,X 3 , ., x n ) be a differentiable 

Lyapunov function such that, 

(a) V{t, x p x 2 , x 3 ,., x n ) > £-(xj,x 2 , x 3 ,., x n ) > 0, 

e (0,0,0,., 0) = 0, V (t, 0,0,0,., 0) = 0. (6.8) 

That is V has a strict minimum at the origin and £ is a 
continuous positive number. 

(b) Now if 

r." 4-+Z--VU <6.-) 

dt [St ; dx, j 
for t > t 0 , then the unperturbed motion is stable. 

3. Outside an arbitrarily small neighborhood of the origin if 

Ytf ^P 2 >0, t>T o >t 0 (6.10) 

and in addition to equations (6.7) and (6.9) if we have 

dV 

-<-/?<0, (6.11) 

dt 

where p~ and [3 are the positive constants, then 


X,. =0(1= 1,2,3, —, n ) is asymptotically stable. 

4. Further, if 

dV 

V>0, V(t, 0,0,0, • • ••, 0) s 0, — < -0 < 0 (6.12) 

dt 

and the derivatives dV / dx j are all finite, then the 
unperturbed motion is stable under constantly acting small 
disturbances. This is Malkina Theorem [6, 8, 21], 


B. Generalized CTTSIP Setup 

The gist of the generalized setup of CTTSIP is 
described below[l-4]. 

Let Y® (t) a be the thermodynamic coordinates of the given 
process on the unperturbed trajectory and K(() a be that on 


the perturbed one. The perturbation coordinates, CC i (t ) a are 
defined as, 

a i (t) = (Y i (t)-Y?(t))^0, (i = 1,2,-", n) (y 

and hence the equation of unperturbed trajectory is 
obtained as. 


a°= 0. 


(7.2) 


Thus on using GC^t) the thermodynamic perturbation 


space gets defined. The constitutive equations of motion 
within the perturbation space in general are described as, 

f.(a l ,a 2 ,a 3 ,---;t) (i = 1,2, •••,«). 

df (7.3) 

The domain of thermodynamic perturbation space is 
determined by 

t>t 0 , t> 0, a t <£, £>0 


and the trivial solution of equations (7.3) is (7.2). 

In CTTSIP[4] the thermodynamic Lyapunov function, JZ S , is 
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defined as, 

4=(<r,(f)-o?(0)^0, (74) 

where, CT (t) and a s (?) are the respective entropy 
source strengths and are positive definite quantities [12, 25, 
26] as per the second law of thermodynamics that appears in 
the Clausius-Duhem inequality [12, 22-24]. Since we are 
using the thermodynamic function, namely the entropy 
source strength in defining the Lyapunov function the 
CTTSIP inherits thermodynamic base. Notice that, C s is the 
excess rate of entropy production per unit volume. 
Obviously, by definition we have, 

4 = C s (a l ,a 2 ,a i ,---,t)>0 


terms in the expression of tr (t) and C s (t) can be 
ignored and hence the expressions for them get 
restricted to. 



a -r( dC ’\ 

{daj 

1 


4 = L s {a x ,a 2 ,a^---,t) < 0, 

and on the unperturbed trajectory we have, 

4=4(o,o,o,-,o = o, 


(7.5) 

(7.6) 

(7.7) 


where C s is a differentiable function [28], As we need to 
work within the thermodynamic perturbation space let us 
expand <T S (t) in terms of CC i (t) a that gives. 


■fe.] 

0 f , \0 

1 v" f d~a s ] 

c?4 

_dC s 

‘{da,} 

a. h 7 a.a ; 

1 2 xjydapaj J ' J 

dt 

fts 

b 


+-Y 

6 4-? 


d 3 a 


6 f~‘ k I da^ajda. 

Notice that, for the sake of brevity the time dependence of 
quantities has not been shown. The thermodynamic 
Lyapunov function, C, , is then given as (c. f. equation(7.4)), 

^ f da ^ 1 f d 2 a s ) 

i ^ da t J Sa t daj J 

1 y- f 3 3 f7 'j 

+ — > - £ - a,a ,a. +.>0. 

6 4-? \ An An, Fir, \ ' 1 k ' 




(7.12) 

In equation (7.12) each (dC s / da t ) may be either 
positive or negative. The preceding assertion gets elaborated 
on differentiating equation (7.12) as follows. 

dC s J d L/gcrV 

8aj ~\da : J “ da, | { da, J 

da s 
daj 

The total time derivative of Jl. then reads ai 


(7.15) 
and hence 
go for the second 


which transforms 

dC s 
dt 


, 


f 5<T 0 


da, 


/ U«J 

Yu' 


if (da s I da j = 0 
{dC s / da, )° = 0, we have 

0 2 er 


order terms, r 


'■ ] *V ' J k J (7.9) 

Alternatively, C s can be expanded directly in terms of a, a 


4 = (a s 


2 ■”[ da i da i 


l^dajdaj) a ’ a ^ 0. 


(7.16) 




da t 


provided (3 2 a s / da.dUj ) ^ 0 . The differentiation of 
4 successively with respect to tt k and a t produces the 
following two expressions, namely: 


biy 

£ 4-! 


da-dajda. 


a i a j a k + .^ 0. 


If (8a s / da, )° i=- 0 and hence 
(dC s / da, )° 0 , and as per the requirement of 

LyapunovQ theory of stability of motion [5-8, 21] 
a, & need to be sufficiently small the higher order 


1 d f \ ’ (da ] 

=-> -— a,a, 

2 da, da,da,. 1 ' J 


i\da,da k 


vf 

' d 2 C 

da t da k 1 

K dafia k y 


ydafia. 
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The total time rate of then reads as 

dC 1 ^ 

—- = > a,a ■ 

dt 2 ff 8t\ dad a, 




r 

ccj 


3. However, if {dcr s I d(X^) =0 


(d 2 cr s !8a t ,da^ = 0, 


we have to go for a third order approximation, 
namely: 


£ s =(cr s -cr s °)~-^X 


f d 3 cr 


6 ~~t l 8a i da ] 8a k j 

(7.20) 

provided ^3 3 cr, / dafia -8a k j ^ 0. The total time 
derivative of £ s then reads as, 

dC 1 8 f 8 3 (T 

—- = — / oc.aM, — -— 

dt 6 ~~^ k 8ty8a i 8a ] 8a k 

-1 f d 3 cr 


~‘\dafia j 8a k 


4. In this sequence if up to (n-l) th derivatives of <J S 
happen to vanish at the origin, one needs to go for 
the n 111 derivative of (T s at the origin. That is, if, 

(8cr s / dafi 0 = 0 , [d 2 cr s I Safi a j ) = 0 and 
so one (( 8 n a s / 8afia j 8a k ■ ■ ■ 8a (n _ ]} )° = 0), 
we have to go for n order derivative, namely: 

C s - (er — <t°) — 

1 ^ f 3V s T (7.22) 

n ! ijt-y8afia j 8a k ■■■ J 

provided (d"cr s /8afiajda k •••j ^0 . The total time 
derivative then read as. 


dC s 
dt n\ 


1 . 8 f 8 n cr 

^ 7 dt\dafiafla h > 


^21 


\dafiajda k ■ 


(a i a j a k •••) ■ 


The conclusions about the thermodynamic stability 
of a process is then drawn using Lyapunov and Malkina 
theorems [5-8, 21] by establishing the signs of £, s , dC, s / dt 


and whether the magnitude of each gradient (8C S / dafi is 


finite or not. 
1. The 


ipcrturbcd trajectory is said to be a 
thermodynamically stable one when, 

c/-±< o. 

dt (7.23) 

The thermodynamic asymptotic stability is obtained 
if, either 


C s >0, 

^<-A< o 



dt 

(7.24) 

C, <0, 

^>A>0, 



dt 

(7.25) 


where /3 is strictly a positive number that vanishes 


only at the origin. 

3. The thermodynamic stability under constantly acting 
small disturbances, as per Malkina theorem, is 
obtained if in addition to equation (7.24) or (7.25) 
each (8fi / dafi is finite. 
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